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AMLra£l 

This paper describes and illustrates two ways of 
performing time-correlated gust-load calculations. The first is 
based on Matched Filter Theory; the second on Random Process 
Theory Both approaches yield theoretically identical results and 
represent novel applications of the theories, are computationally 
fast and mav he applied to other dynamic-response problems. 
A theoretical development and example calculations using both 
Matched Filter Theory and Random Process Theory approaches 
are presented. 


Nomenclature 

Roman 

G(s) transfer function of vertical velocity component 
of von Kannan turbulence (gust) 

hi(t) time response of output variable i to unit im- 
pulse 

I Vi(co) frequency response function of output variable i 
K arbitrary constant 


w g (t) vertical velocity component of gust 

W g (o>) Fourier transform of w g (t) 

x(t) "matched" excitation waveform 

X(0)> Fourier transform of x(t) 

y(t) time response of output quantity y 

y (t) time response of output quantity y to excitation 

' * matched to y 

y 7 (t) time response of output quantity y to excitation 
matched to z 

Y(o>) Fourier transform of y(t) 

z(t) time response of output quantity z 

z y (t) time response of output quantity z to excitation 

matched to y 

z z (t) time response of output quantity z to excitation 
matched to z 


L scale of turbulence in G(s), equation ( 1 1) 

Rj (t) cross-correlation function between quantities ) 
' and j; auto-correlation function if i equals j 

5 Laplace variable 


Z(co) Fourier transform of z(t) 

Greek 

Cj rool-mean-square of quantity i 


Sij(co) cross-power spectral density function between 
1 quantities i and j; auto-power spectral density 
function if i equals j 

t time 


time argument for auto- and cross-correlation 
functions 

circular frequency (rad/sec) 


t 0 arbitrary time shift 

V velocity of aircraft in equation (1 1) 


Superscripts 
complex conjugate 
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Subscripts 


g gust 

h unit impulse time response 

x matched excitation waveform 
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y response quantity y 

z response quantity z 


Introduction 

The primary structure of an aircraft must withstand all 
the static and dynamic loads it is expected to encounter during its 
life. Such load's include landing and taxi loads, maneuver loads, 
and gust loads. There are several methods used by airframe 
manufacturers to compute gust loads to satisfy certification 
requirements and include both stochastic methods and 
deterministic methods. The US. Federal Aviation 
Administration (FA A) recently asked the National Aeronautics 
and Space Administration (NASA) for assistance in evaluating 
the Statistical Discrete Gust (SDG) Method [1| as a candidate 
gust-loads analysis method for complying with FAA certification 
requirements. The SDG method is a time-domain approach and 
yields time-correlated gust loads by employing a search 
procedure. 

During the course of the NASA evaluation of the SDG 
method* the authors recognized that Matched Filter Theory 
(MFT) [21 could be applied to the linear gust problem to 
compute time-correlated gust loads. Computing time -correlated 
gust loads in this manner has the twin advantages of being 
computationally fast and of solving the problem directly. 
Historically* MFT was first utilized in the detection of returning 
radar signals, Papoulis (3) showed that the principles of MFT 
can be used to obtain maximized responses without the use of 
calculus of variations and in applications other than signal 
detection. The first purpose of this paper is to demonstrate that 
MFT is also applicable to aeroelastic systems, and specifically 
applicable to the computation of time-correlated gust loads. 

During the course of the MFT investigation, the authors 
found that the maximized and as well as the other responses of 
the system looked similar to auto- and cross-correlation 
functions of Random Process Theory (RPT) [4|. From this 
observation and also from theoretical considerations given later* 
it was recognized that time-correlated gust loads computed by 
MFT are theoretically identical to those obtained using RPT. To 
the knowledge of the authors* the correlation functions of RPT 
have* until now, not been interpreted as time-correlated gust 
loads. Some preliminary analyses demonstrating this 
interpretation were presented by the authors in reference 5. The 
second purpose of this paper is to show the relationship between 
the results of MFT and RPT. 

Both the MFT and the RPT procedures for computing 
time-correlated gust loads involve novel applications of the 
theories and unconventional interpretations of the intermediate 
and final results. This paper outlines the mathematical 
developments* recognizes a duality between MFT and RPT* and 
presents example calculations using both MFT and RPT for 
computing time-correlated gust loads. 


Backsrgmid 

This section of the paper defines time-correlated gust 
loads and the current analysis capabilities in this area. 

Time-Correlated Loads 

This paper deals specifically with time-correlated gust 
loads and figure \ illustrates two types of such loads. Time- 
correlated loads are time histories of two or more different load 
responses to the same disturbance quantity. In the figure the 
disturbance quantity is the vertical component of one- 
dimensional atmospheric turbulence and the time-correlated 


loads (the output quantities) are the resulting load responses (in 
this case the bending moments and torsion moments) at several 
locations on the airplane wing. 

The first type of time-correlated load is illustrated on the 
right wing: loads (two bending moments in this illustration) at 
two different locations on the airplane. The second type is 
illustrated on the left wing: two different loads (bending 

moment and torsion moment in this illustration) at the same 
location on the airplane. 

As indicated by the time histories in the figure, time 
correlation provides knowledge of the value (magnitude and 
sign) of one load when another is maximum (positive or 
negative). In general, the time correlation principle allows the 
computations of values of all loads at all the stations with respect 
to the time when any one maximum load occurred. Such 
information may be used directly during analyses and testing of 
aircraft structures [6]. 
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Figure 1* Time-Correlated Gust Loads 


^S tatistical Discrete Gust Methud 

The SDG method determines the response time histories 
of "worst case” gust loads (such as shear forces* bending 
moments, and torsion moments) and the corresponding critical 
gust profiles" which produce them. These loads are lime 
correlated and this feature is a major advantage of the buo 
method over some other gust loads analysis methods. 

Another advantage of the SDG method is its applicability 
to nonlinear* as well as to linear, systems. The SDG method 
employs a search procedure to obtain its answers. For a linear 
system, by taking advantage of superposition* the search 
procedure may be simplified and reduced to a one-dimensional 
search For a nonlinear system, however* this is not the case 
and the resulting search procedure remains multi-dimensional 
and can become computationally expensive [7]. 


Scop e 

In light of the technical contributions provided by the 
SDG method as outlined above* this section states the conditions 
under which and the means by which this paper makes a 
contribution to the area of time-correlated gust-load calculations: 
(a) the SDG method is capable of performing both linear and 
nonlinear analyses* whereas the present methods are restricted to 
linear systems only; (b) the SDG method uses a search 
procedure, whereas the present methods obtain the same 
quantities directly and with a significant reduction in 
computation time. 
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T>poreticnl Develop ment 

In this section the theoretical background of the MFT and 
RPT approaches will be presented. For this paper, the Fourier 
transform of a function f(t) is defined by, 

F(U)) = J f(t) e (it, 

and the inverse Fourier transform is given by, 

f(,) = 2irj F(C0)eiuM d0) - 


Maximum Respons e tn "Matched" Excitation Waveform - 
Equations (la) and (lb) express the fundamental result of 
Matched Filter Theory and figure 2 illustrates the result: the 
response of the system. "filter" to a unit impulse is proportional 
to the desired input waveform, x(t) f when delayed by t 0 and 
reversed in time. For convenience and to be consistent with 
references 3 and 8, the form given by equation (la) will be used 
in the remainder of the theoretical development of the MFT 
approach. It should be pointed out that either form of the the 
equation will lead to the same final results. 

The system response to x(t) may be obtained from the 
inverse Fourier transform 


V Approach 


The objective of MFT, as originally developed, is the 
desicn of nn electronic filter such that its response to a known 
input signal is maximum at a specific time [3,8|. It found early 
application to radar considering the 'Tiller" to be a detector that, 
in response to a known input signal, produces an output signal 
for further processing [2|. In this case the detector design is the 
optimum design for maximizing the output signal-to-noise ratio. 
The matched filter is designed so that the unit impulse response, 
h (t), of the output y of the filter is proportional to the known 
signal, shifted and reversed in time, 

h y (t) = Kx(-t + t 0 ), Oa) 

where K is an arbitrary constant, t 0 is an arbitrary time shift, and 
x(t) is the known input signal. As will be shown later, the time 
shift may be selected through criteria of convenience only and 
will be the time that the matched filter response is maximum. 

In the present application, the filter is considered to be a 
system whose dynamics are known. Specifically, the system is 
characterized by the combination, in series, of the dynamics of 
atmospheric turbulence and the dynamics of aircraft load 
response. Thus h y (t) is known. The input signal. x(t), is lo be 
determined from the impulse response of the system. From 
equation (la), 

x(t)=hy(t 0 -t)/K. (lb) 


y (t) = JL J Y(w) e dw, (2a) 

where 

Y(w) = H y (co) X(oa). (2b) 

and H y (o>) and X(w) are the Fourier transforms of h y (t) and x(t)- 
H y (co), which is the frequency response function of the output 
y, may be determined from equation ( la) as 

H y (o>) = Ke X(-m) 

= KX*(co) e i0 *" . (3) 

Thus the response y(t) is 

y(t) = J- J KX*(o>)X(a>) e "° <l dco (4) 

•«* 

Note that the product X*(co) X(co) is the energy spectrum or 
power spectral density (PSD), S„( w )i °f *(0. and is an even 
function of frequency, so 

y ( t) = K ^ J S^fcoXosUiXt - 1 0 )) dw ( 5a ) 


The simple result of MFT allows direct determination of the 
input signal, or excitation, that produces a maximum response of 
the system. The result, as will be shown, is the maximum load 
response and the critical gust profile. 


The maximum response occurs when t - to, 

y m « = y('o> = K[^Js„(tt)d(o] 


(5b) 


Unit Impulse Impulse Response 





Matched Filter Theory 
Result 


r lull = Kx(* M - 1„)> 
Mill t„-n 




Matched Excitation Waveform 


x <t) 



H F 


-► y (*) 


Response to Matched 
Excitation Waveform 



Figure 2. Matched Filter Theory - Signal Flow Diagram 


It is convenient to set the term in brackets in equation (5b) equal 
to unity, which means that the root- mean- square (RMS) of the 
excitation waveform, x(t), is unity. That is, 

-L f x‘((0) X(to) dot = J- f S„(co)dco= I (6) 
2k J 2*2 

It can be shown that a necessary result of this assumption is that 
the arbitrary constant, K, is equal to the RMS of the impulse 
response, a», y • The stipulation in equation (6) that the 
waveform be of unit energy will be used later as a constraint on 
all excitation waveforms considered and will provide a basis for 
comparison of their effects. 


Finally, substituting for K in equations (5a) and (5b) 

gives. 
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yU) = a h 2 ^ J s x*( w ) cos(0)(t - t 0 )) da) (7a) 

and 

y,nax = XU = \ * 7b) 


and H y (0)) is the frequency response function of an aeroelastic 
system to turbulence. With the turbulence ’’pre-filter” known, 
the critical gust profile can be obtained as an intermediate result 
as 

W g (co) = G(co)X(cd) 
and from equations (3) and (13), 


The input x(t) is said to be "matched" to y, since x(t) is 
determined from the impulse response of y , equation (1), and 
will be referred to as the matched excitation waveform. 


Response to Arbitrary Excitation Waveform - It can be shown 
that the response y'(t), say, of y to any arbitrary excitation, say 
x'(t), that satisfies the normalization constraint of equation (6), 
will never exceed the maximum, equation (7b), produced by the 
matched excitation, x(t). To see this, the Fourier transform of 
this response may be written 

Y’(co) = H y (co) X'(o>) . (8) 


Substituting equation (3) for H y (co) and taking the inverse 
Fourier transform gives 


y’(t) = | o hy X*((o)X f (co) e ^ tf>) dto . 


(9) 


Applying Schwartz's Inequality to equation (9) gives 


|y'(t)| <> cl 


-L f X*(co)X(co) da) 
2 k J 

| X , *(w)X’(co) da) 


X 

2 n 


(10a) 


Both bracketed terms equal unity because of the constraint, 
equation (6), imposed on the waveforms. Thus 


W g (w) = G(co)G’(a))H‘(co)e ( 14a) 

h y 

and taking the inverse Fourier Transform, 

W » (0 “ij w * ( “ )ci " d “ 


G(o))G*(a))H y (o))e 1<0<l Io) da) (14b) 


This theoretical development was carried out in the 
context of the response of an aircraft to atmospheric turbulence, 
where the pre-filter represented the atmospheric turbulence. 
Theoretically, other pre-filters can be used to analyze responses 
to some other disturbance quantity. Potential alternative pre- 
filters include runway roughness, pilot dynamics, and control- 
surface deflections. The range of alternatives is limited only by 
the ingenuity of the analyst in developing appropriate pre- filters. 



Duality between MFT and RPT 

Consider now two excitation waveforms x y (t) and x 7 (t) 
that are matched to two outputs, y and z, of a given system. The 
response of y to x z (t) is then 

y 7 (t) = °h, } Xy(<o)X z (o)) c ,a><l ' o) dca . (15) 

From equation (3), 


|y'(t)U<r . (iob) 

hy 

which states that the magnitude of the response of y to any 
excitation waveform (suitably normalized to satisfy equation (6)) 
never exceeds the maximum response to the excitation waveform 
that is matched to y. 


Any given frequency response function 

can be built up from several separate ones and thus represent the 
combined effects of several subsystems. In the context of the 
present paper, a von Karman turbulence "pre-filter”, represented 
by a transfer function, is implemented ahead of the aircraft 
response function. The turbulence transfer function used in this 
study is adapted from that obtained from reference 9, 


G(s) - O, 


(1 + 2.618 *) (I * 0.1298 

(1 + 2.083 s) (1 ► 0.823 ^ s) <1 4- 0.0898 s> 


(in 


Thus, the transform of the response y, as expressed above in 
terms equation (2b), would be rewritten as, 

Y(o» = H y (to)G(o)) X(o)) (12) 

where 

H y (co)G(co) = H y (a)) (13) 


and 


X y ((D) = 


H;(co) . 


i<ot 0 


Thus, 



- icx* 

e 


(16a) 


(16b) 




(0)) e 


Mi u 


da) . 


(17) 


The product H y (co) H z *(a)) is the cross-PSD function. S yi (o>), 
for outputs y and z, thus. 




S y *(o>) c 


Mi - U 



(18) 


The term in the brackets is the cross-correlation function, R yz (t - 
U, between outputs y and z [10), so. 


y,(0 = — R yI (t -><,)■ 

h t 


(l 9a) 
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Similarly, the response of z to x y (t) would be 

z y (t)=J-R, y (t-t„), 09b) 

Hy 

where 

if * i«K* - l o) . 

R /y (t - t n ) = -^ I M,(co)H y (to) e dco 

= J_ fs 7y (( 0 )e i “" U dfi>. (20) 

2 It J ^ 

Note also that these two cross-correlation functions satisfy the 
relation (where t = t - to )» 

R ij (t) = R ji (-t). 

Notice that if the response of y (or z) to its own matched 
excitation, x y (t) (or x 7 .(t) ) were being considered, auto - PSDs 
(S vv or S/.,.) and auto - correlation functions (Ryy or R zz ) are 
involved. Thus, dynamic responses to various matched 
excitations are proportional to auto- and cross-correlation 
functions between the outputs. Further, when t = t 0 , 


y y (' 0 ) = a~ R y/ 0) ' (21) 

N 

The auto-correlation function is an even function of its time 
argument, x, and its value is maximum at x » 0 and is equal to 
the mean square of the PSD, Syy. Thus, 

y y (t 0 ) = y m .* = \ • (22) 

which is Ihc result obtained from the MFT approach, equation 
(7b). 

The duality between MFT and RPT is apparent from 
comparison of equation (15), which results from the application 
of MFT with equation (19), and from comparison of equation 
(7b) with equation (22). These equations show that the 
correlation functions of RPT can be interpreted as scaled time 
responses to the "matched" excitation waveforms of MFT; and, 
necessarily, that the value of the autocorrelation function at X - 
0, R (0), is proportional to the maximum response, y m ><> of 
MFT The duality between MFT and RPT also allows a 
relationship to be established [11] between MFT and Phased 
Design Loads Analysis [61. Further, since the product 
G(o))G*(w) in equation (14b) is equivalent to the von Karman 
turbulence PSD, the critical gust profile resulting from the 
matched excitation waveform can also be determined from 

w g (t) = — — - — [ <D g (co) H y (co)e ,<0(l '° ) dco . (23) 

n y .i 

Thus, RPT, normally employed to determine the statistical 
properties of random processes, may also be employed to obtain 
certain deterministic responses provided that they axe interpreted 
as responses to matched excitation waveforms as described by 
MFT. In circumstances for which no PSD is available (e.g. 
maneuvers) the MFT approach would have to be used. Actual 
implementation of the two approaches will be discussed 
presently. 


Aj mlkalmu or the Analysis Juimliirgs 
lo Time-Cor rel ated (iust Loads 


MFT Approach 

Figure 3 contains a signal flow diagram of the steps 
necessary to generate a maximum dynamic response at some 
point in the aircraft structure. It is the mechanization of the 
information presented in the Theoretical Development section of 
the paper. The signal flow diagram is presented in two parts; 
the top half illustrates the generation of the system impulse 
response; the bottom half illustrates the generation of the 
maximum response of the system. 

In the top half, the gust spectrum is excited by an 
impulse of unit strength to generate an intermediate gust impulse 
response which, in turn, is the excitation to the aircraft. 
Computationally, the time history of the impulse response is 
carried out until its magnitude dies out to a small fraction of the 
largest amplitude. Several output channels (y, z\, ...) arc 
shown in the Figure, however, only one response (y) is used. 
This impulse response is then normalized according to equation 
(6a). 

The bottom half illustrates how all the time-correlated 
responses are obtained including the maximum response of the 
output y, and the critical "gust" profile. The y response builds 
to a maximum at t = t 0 at which point the excitation ends. The 
response then decays to near zero. The maximum response, 
y is equal to the RMS of the system impulse response. It 
should be mentioned that the responses shown for the other 
outputs are not maximized. The critical gust profile and 
maximum response of the system are unique to only one load 
output (the top output, y). For other maximum load responses a 
separate but similar analysis needs to be performed. 

An important detail illustrated in this figure is the 
introduction of the a pre-niter. The effect of the pre-filter is to 
provide dynamics of the input disturbance which itself 
contributes lo the shape and magnitude of "matched excitation 
waveform. In this example, the pre-filter is an s-plane 
approximation of the von Karman spectrum, but in other 
applications it could be landing or taxiing disturbance dynamics 
or possibly "pilot" dynamics for maneuver loads. 


RPT A do roach 

Figure 4 contains a signal flow diagram with two parts 
and is analogous to the diagram in figure 3. This figure 
illustrates the steps necessary to generate time-correlated gust 
loads using RPT. Whereas the signals in figure 3 were all in the 
time domain, all but one of the signals in this figure are in the 
frequency domain. 

From the top half, the "Known Dynamics" box is the 
same as that in Figure 3. The input to this box is Gaussian white 
noise and the output includes auto- and the cross-power spectral 
density functions of some aircraft response, with an intermediate 
output being the von Karman power spectral density function of 
atmospheric turbulence. 

Time-correlated gust loads arc obtained in the bottom 
half of the figure by taking the inverse Fourier transforms of the 
auto- and cross-power spectral density functions obtained in the 
top half. It should be mentioned that to obtain precise 
representations of the time-correlated loads it was necessary to 
deal numerically with two-sided spectra (that is, with both the 
positive and negative frequency components present). 
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Implementation of Matched Filter Theory Approach - Signal Flow Diagram 
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Figure 4. Implementation of Random Process Theory Approach - Signal Flow Diagram 








IMnineric fll Examples 


f \prnpinsiic Configu ration 


S? 

II - - \v- 


In the application of MFT and RPT, existing structural 
and aerodynamic models of the NASA DAST A R W' 2 were 
used This configuration is a Firebee II target drojie fitted w 
an Aeroelastic Research Wing (ARW) [I2|. This model was 
especially well-suited for the study since it has structura 
flexibility, a stable and dominant short period, and several load 
outputs. The load outputs are comprised of shear forces, torsion 
moments, and bending moments at several points along the span 
of the wing. A detailed description of the dynamic loads model 
is available from reference 13. 


Analysis Conditions 


Altitude = ! 5,000 
Mach No = 0.7 


Vehicle Description 

Wciqbl = 2500 1b 
Winn r>pnn - 227.04 in 
neference Lenqlh = 2'1'1/in 
Aspect n,it»o = 10 3 


Annlyhc.il Representation of Vehicle 

2 Lonqittntinal HiqiH-nody Modes 
U Symmeliic flexible Modes 


Analytical Representation ol Turbulence 

von Kami an Speclrum 
I. = 2500 H 
= I in/soc 

W <1 


Figure 5 shows some of the more important data of the 
vehicle itself and of the analytical representation of atm ° s P^ 
turbulence. Tlic structural part of the model was derived from i a 
finite element code and the unsteady aerodynamics (at a Mach 
number of 0.7) from a doublet lattice code. Besides the two 
rigid-body modes, eight symmetric flexible modes ^ ere ^' aine ^ 
for this study. The final dynamics equations (the quadruple 
equations), constructed with a matrix analysis code, consisted o 
97 first order equations, 9 output equations and one inpu . 
These final equations contained the dynamics of the structure, 
the of the unsteady aerodynamics, the loads, and the von 
Karman spectrum. 


Time-Correlated GusLio nri s Using Matched Filia l 
Theory 

Fieure 6 gives shows of the time-correlated loads at the 
wing root and the wing tip stations. The top two plots on the 


Figure 5. Aeroelastic Configuration - NASA DAST ARYV-2 

righthand side contain the time -hi stories of bending moment and 
torsion moment at the same root station on the structure. The top 
plot contains the bending moment response resulting from the 
excitation matched to root bending moment. The critical gust 
profile matched to the root bending moment is contained in the 
left plot. Although their magnitudes and signs may be different, 
the general shape of waveform represented by this plot is very 
similar for almost all the critical gust profiles of the system. The 
critical gust profile responses have the characteristic of being 
dominated by the short period dynamics of the aircraft. This 
phenomenon can be explained by examining equation (14a). 
The higher frequency aeroelastic responses due to the elastic 
modes are doubly attenuated by the high frequency roll-off of 
effectively two von Karman ’’pre-filters.” The bottom righthand 
plot shows the time-history of bending moment response at a 
location near the tip of the wing. 



Figure 6. Time-Correlated Gust Loads - Matched Filter Theory Approach 
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Root Bending Moment 


Root Torsion Moment 




Time, sec 


Time, sec 


Figure 7. Comparison of MFT and RFT Calculations 


Tomnarison of Res u lts From MFT and R . E I 
Calculations 

Figure 7 shows a comparison of wing root bending 
moment and torsion moment time responses calculated with the 
MFT and RPT approaches. The Matched Filter Theory results 
are the same as those shown in figure 6. For purposes of 
comparison, the s-plane approximation of the von barman gust 
power spectral density function is used for both the MFT and 
RPT calculations. Except for some slight differences in the 
peaks and troughs, results from the two approaches are 
practically indistinguishable. This is in accordance with the 
MFr/RPT duality described in the Theoretical Development 
seciion of the paper. 

This paper has described and illustrated two approaches 
for computing time-correlated gust loads. The two approaches 
vield theoretically identical results and are computationally fast, 
the first is based on Matched Filter Theory and is a time-domain 
approach; the second is based on Random Process Theory and is 
a frequency-domain approach. These approaches represent 
novel applications of the theories and unconventional 
interpretations of the intermediate and final results. 

An important contribution of this paper is 'he 
introduction of a "pre-filter" in the application of Matched Filter 
Theory. The pre-filter is the mechanism by which the dynamics 
of the input disturbance (atmospheric turbulence) are 
incorporated into the problem and the resulting matched 
excitation waveform contains these dynamics. The output of the 
pre-filter due to the matched excitation waveform is the critical 
gust profile. 

A second contribution of this paper is the interpretation 
of scaled correlation functions from Random Process Theory as 
time histories. It has been shown that the time-corTelated gust 
loads predicted by Matched Filter Theory are theoretically 
identical to the auto- and cross-correlation functions, 
appropriately scaled, predicted by Random Process Theory 
Further, the critical gust profile obtainable from the Matched 
Filter Theory approach is also obtainable from the Random 
Process Theory approach. 

Although, in this paper. Matched Filter and Random 
Process Theories have been applied to the calculation of tune- 
correlated gust loads, the approaches are general enough to he 
applied a variety of dynamic-response problems, such as taxi 
and landing loads, maneuver loads, and gust loads. The choice 
of which to use depends on the intended application. 
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